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The perturbative QCD calculation of heavy quark fragmentation into a heavy 
meson is extended to predict fragmentation into heavy flavor baryons. This is 
accomplished by implementing the quark-diquark model of the baryons. Several 
diquark form factors are used to enable the integration over the virtual heavy 
quark momentum. The resulting spin independent functions for charmed quarks 
to fragment into charmed baryons with spin 1/2 and 3/2 are compared with recent 
data. Predictions are made for the spin dependent fragmentation functions as well, 
particularly for the functions g\ in the case of spin 1/2 baryons. 



1 Introduction 

Fragmentation functions have received considerable attention in recent years. 
While experimental information beyond the pion distribution (presumably 
from light quarks) has been slow in accumulating, theoretical interest has been 
growing. The particular functional form for heavy flavored quarks to frag- 
ment into heavy mesons has been studied using Operator Product Expansion 
techniques, light cone quantization, QCD perturbation theory. Heavy Quark 
Effective Theory, and other methods. Some of these methods yield general 
properties that reflect the overall structure of QCD, as it is currently under- 
stood. Other approaches take particular models of the low energy behavior 
expected from QCD, but in regions that are not perturbatively calculable. 
The particulars of the various approaches will be subject to some experimen- 
tal scrutiny in the future, but are not yet put to the test. One general feature 
is known - the peak of the hadron distribution moves toward higher momenta 
as the quark mass increases. This feature is a result of the kinematics im- 
plicit in most models of the non-perturbative process, and is incorporated in 
the phenomcnological Peterson functiorEl that is used to fit the sparse data on 
heavy quark fragmentatiom. 

Even more difficult to test experimentally are the spin dependences of the 
fragmentation processes. Yet these dependences are important to know. They 
reflect the details of the primarily non-perturbative mechanism by which par- 
ton polarization is passed on to the hadrons. In this sense, the spin-dependent 
fragmentation involves the reverse of the process by which the nucleon spin 
is shared by its partons (the "spin crisis"), and may reveal a similarly myste- 



1 



rious decoupling of valence quark spin and liadron spin for some regions of 
kinematics. 

Over the last several years a number of theorists have noticed that for 
fragmentation of a heavy flavor quark into "doubly heavy" mesons, like the 
c-quark into the J/vt or the b-quark into the Be, perturbative QCD may be 
applicablflu. If this is the case, the fragmentation functions are calculable, 
at an appropriate scale, and QCD radiative corrections can be obtained from 
the Renormalization Group or the Altarelli-Parisi equations. Such calcula- 
tions have been performed and scrutinized. It has been shown that in the 
heavy quark limit (i.e. the mass goes to infinity) the functions have the form 
expected from more general consideration^. This corresponds to the heavy 
meson taking all of the heavy quark's momentum; the distribution becomes a 
delta function at z = 1. The l/mq corrections are calculated alsoQ. In any 
case, this approach can predict the spin-dependent fragmentation functions 
along with their momentum and mass dependences. In the heavy mass limit, 
of course, the spin of the heavy quark is conserved, so the spin dependence is 
simple. What is of phenomenological interest is the next order correction, at 
least, since that has non-trivial spin dependence. 

In some circumstances the spin dependence of fragmentation is most read- 
ily studied experimentally by observing baryons rather than mesons. This is 
true for the production of hyperons or heavy hyperons (Ac, Ab^tc), wherein 
the weak, parity violating decays provide polarization analysed. To consider 
fragmentation into baryons in this perturbative scheme, the three quark sys- 
tem has to be confronted. A simple alternative is to consider the baryons as 
quark-diquark bound states, and to use the same perturbative method as for 
the mesons. In order for the perturbative calculation to be useful the creation 
of a heavy pair of quarks or diquarks must be an intermediate step. Hence, 
doubly heavy baryon fragmentation is an appropriate testing ground for these 
ideas. It is not expected that sufficient data to study this process will be 
available in the near future, however. 

To begin to see the structure it will be worthwhile to stretch the region 
of applicability to the "singly" heavy baryons. We have been carrying out 
this program to see the expected spin and kinematic dependences, with the 
hope of providing an experimentally testable model. Of immediate interest 
is the question of whether the baryon fragmentation functions have the same 
kinematic dependence as the meson case. In general the answer is no in this 
model, but the details have to be studied. Furthermore, the spin dependent 



"^It is not clear to the present authors that all sources for such terms have been consid- 
ered, particularly corrections coming from the treatment of the Bethe-Salpeter bound state 
wavefunction for the meson 
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fragmentation is interestingly distinct from the naive heavy quark limit. The 
perturbative calculation and its results will be presented below, along with a 
comparison with some recent data. 



2 Perturbative calculation 



The first calculations of the fragmentation functions in the perturbative scheme 
were applied to some of the inclusive heavy flavor meson decays of the as 
produced at LEP. The partial width for this inclusive process can be written 
in general for a hadron H as 

dr{Z° ^ H{E)+X) = y^ [ dzdllZ'^ ^i{E/z)+X,fi)Di^H{z,ii), (1) 

. Jo 

where H is the hadron of energy E and longitudinal momentum fraction z 

relative to the parton i, while fi is the arbitrary scale whose value will be 
chosen to avoid large logarithms. The fragmentation function enters here in a 
factorized form (that can be maintained through the evolution equations). 

Now, consider the final state with one heavy flavor meson, say the Be for 
definiteness. To leading order the Be meson arises from the production of a 
pair of b-quarks, in which one of the quarks fragments into the meson. As 
Fig.l illustrates, with Q = b,Q' = c,Q' = c, the perturbative contribution 
involves the virtual b- 
quark radiating a hard 
gluon (in axial gauge, 
so that there is no con- 
tribution from the op- 
posite quark) . The 
hard gluon produces a 
heavy flavor pair of c- 
quarks. The gluon 
must have energy at 
least twice the charm 
mass, so that the cou- 
pling asiQ"^) is small. 
For matching momenta 
(or relative 3-vclocity 

zero) the b and c form the Be meson with probability given by the square of 
the Bethe-Salpeter wave function. Since the doubly heavy mesons are weakly 
bound objects, the wave function at the origin is known from non-relativistic 
quark models for the heavy-heavy system. 




Figure 1: The amplitude for Z° 
Baryon(QD)+X. 



MesoniQQ') + X or 
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The amplitude for Fig. I, Ai, can be evaluated explicitly from perturbation 
theory. The decay rate for unpolarizcd Z'^ ^ Be + c + b can be written 

Ti = ^Jm[dp]W]i2nrS*{Z-q-p-p')^J2\'^,\^ (2) 

where q, p, and p' are the 4-momenta of the b, Be, and c, respectively. 

To obtain the full inclusive width the unobserved quarks must be inte- 
grated over. The phase space integration can be simplified considerably when 
the limiting case of Alz — > oo is approached by taking leading order in rrib/Mz 
(with nic < mi,). The two body phase space for p and p' can be written as an 
integration over z and s, with transverse momentum fixed for each such pair. In 
the large AIz approximation the transverse momentum of the hadron is small 
and p = zq. Once the square of the amplitude Ai is summed over spins and 
simplified by dropping non-leading contributions, the width for Z^ — > cc can be 
factored out of the expression Eq. |l] leaving an integral over the fragmentation 
function. Then 

/ dzDc^Hiz) ^ '1 ^ ^' / ds dzQ{s ^---^)F{s,z), 3 

Jo 27'Kmc Jo Jo z 1- z 

where R(0) is the Bethe-Salpeter wavefunction at the origin, and F is the 
remaining integrand, which depends on s — q^, z and the quark masses. The 
upper limit on the s integration appears as the Mz oo limit. So the partial 
width for Z'^ —> H + X is given by an integral over the virtuality of the heavy 
quark and the phase space of the unobserved degrees of freedom. 

The same procedure can be applied directly to the baryons, if the quark- 
diquark model of the baryons is used. The hard gluon in the process must 
produce a diquark-anti-diquark pair, and the diquark (color anti-triplet) com- 
bines with the quark to form the baryon. Note that an alternative scenario has 
the heavy quark fragment into a diquark first, and then the diquark dresses 
itselLto form the baryora. This leads to very different results, as pointed out in 
RcfE3, and is not justifiable hfirein, where the diquark is not necessarily heavy 
fiavor. These latter authors Ii3 have performed a calculation that is similar in 
spirit to part of the spin independent procedure we follow below. 

We now proceed with the calculation of fragmentation functions for (singly) 
heavy flavor baryons. The basic covariant coupling of diquarks to gluons was 
written long agcQ. There is one coupling constant for the scalar diquark color 
octet vector current coupling to the gluon field — a color charge strength, along 
with a possible form factor F^. The momentum space color octet current 
(which couples to the gluon field vector) is 

J^^(^) =5.F,(fc2)(p + p')^^"tAA^^/3^ (4) 
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where p and p' are the scalar diquark 4-momenta and k = p' —p. For the vector 
diquark there are three constants - color charge, anomalous chromomagnetic 
dipole moment n, and chromoelectric quadrupole moment A, along with the 
corresponding form factors, Fe, Fm, andi^g. 

= gs{\'')po^{FE{e)[e^{p)-tP\p')]{p + p')^ 

+ (1 + n)FM{k^)[e'^^{p)p ■ e^Hp') + 6^t(yy . ^"(p)] 

+ ^FQ{k^)K{p)4Kp') + bp.e'^ip) ■ e^Hp'Wk'^ip+p'), } , 

(5) 

where A is the color octet index, a, (3, are color anti-triplet indices, the e's 
are polarization 4- vectors for the diquarks. 

In the perturbative diagrams involved here, the virtual heavy quark emits 
a time-like ofF-shell gluon, that produces a diquark-antidiquark pair, while 
attaining nearly on-shell 4-momentum. The diquark combines with the heavy 
quark to form a heavy flavor baryon, whose amplitude for formation is related 
to the Bethe-Salpeter wavefunction for the diquark-quark system. As in the 
meson production calculations, it is assumed that the constituents are heavy 
enough so that the binding is relatively weak, i.e. the quark and diquark are 
both on-shell and the binding energy is negligibly small. This is expected to be 
true for constituents with masses well above Aqcd, and even the light flavor 
diquarks almost satisfy this constraint. The basic perturbative amplitude is 
shown in Fig.l with the Q'-quark line replaced by an (anti-) diquark D line. 

It should be realized that the integration (over s, the square of the vir- 
tual heavy quark mass) involved in the calculation would diverge for point-like 
vector diquarks, since the gluon coupling to a pair, Eq. ^ carries momentum 
factors. The virtual mass in the integration, ^/s, is passed on to the gluon and, 
subsequently, to the gluon-diquark vertex. Hence it is essential to regulate the 
integrand by some means. This is best accomplished via the chromoelectro- 
magnetic form factors for the gluon coupling to the diquark. The form factor 
approach makes physical sense - it is a result of the compositeness of the 
diquarks. And for consistency, once the vector has form factors, the scalar 
diquark must have one also. 

There is no direct information about the chromoelectromagnetic form fac- 
tors. We may expect that the ordinary electromagnetic form factors will have 
the same functional form as their QCD counterparts — the source of both sets 
of form factors is the matrix element of a conserved vector current operator. 
In the relevent case here, though, the vector operator is the gluon field — a 
color octet. Also, what is of concern here is the time-like region of the form 
factor. For diquarks, of course, there is not any direct empirical evidence 
about their electromagnetic form factors, but diquark-quark models of the nu- 
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cleon have constrained the parameterization of the form factors. Dimensional 
counting rules require that_asymptotically, Fs ^ l/kP, Fe or m ^ Vkl^; ^^'^ 
Fq ~ Jacob, et alEl, have obtained electromagnetic form factors for 

the diquarks (as haa^ recent study of higher twist contributions to the nucleon 
structure function£3). The diquark form factors are assumed to have simple 
pole or dipole forms, and the resulting pole positions appear near 1 GeV. If 
we make the assumption that the color form factors have the same functional 
form as these empirical electromagnetic form factors, we can proceed. In the 
integration that will be done here, the time-hke 9^ region begins at 4m|)j^„^^j, 
(below the NN threshold) for the value z = 1/(1 + rriD/mB), and at higher 
values for other choices of z. This implies that the integration region either 
overlaps or comes near to overlapping the pole positions. We treat the poles 
as real resonance positions, by including a sizeable imaginary part (of 1 GeV). 
This is sensible physically, since the color octet form factor would be dominated 
by color octet vector mesons, and the latter are not expected to be strongly 
bound or narrow. Hence we have the Breit-Wigner forms and their squares, 
with pole positions as given by Jacob, et ai. This choice hides our ignorance 
and provides an interpolation between the space-like and time-like asymptotic 
regions. 

The amplitudes for the baryon production can now be calculated. The 
spin 1/2 ground state baryons are composed of a scalar diquark and a heavy 
quark in an s-state. There is only one coupling, and it involves the Fs- The 
amplitude is 

As 1/2 = -^j^Fs{k^)UBgs[kx - 2mdVx]P\ (6) 



where 



= A^g.r ^^f + 1+^ . (7) 
(s - mA) 



For the vector diquark baryons, there are two form factors (we take the 
quadrupole to be zero - it falls as l/l^l^ asymptotically). The chromomagnetic 
coupling involves a parameter k, the "anomalous chromomagnetic moment" . 
This is taken as 1.39. The s-state baryons are spin 3/2 and 1/2, which we 
will refer to as 1/2'. The 1/2' lies between the 3/2 and the ground state 
1/2 baryon. The amplitude for vector diquarks to be produced, along with 
the heavy quark, contributes to both 3/2 and 1/2' states. The amplitude is 
conveniently divided into a chromoelectric and chromomagnetic part, involving 
the two distinct form factors. The chromoelectric part contributing to the spin 
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1/2' baryon is 



The chromomagnetic contribution to the spin 1/2' baryon is 



(8) 



i-M 1/2 — 



m 



FM(fc2)(l + K)UB75r^9s[9^cx{e*v)md - e*),k^]P^. (9) 



For the spin 3/2 baryon the corresponding amplitudes are 



1E3/2 



m 



and 



M3/2 



^(0) 



(10) 



(11) 



In these amplitudes, V'(O) is the Bethe-Salpeter wavefunction at the origin 
(for the s-state Q-diquark system), is the appropriate diquark mass, Ub 
is a spin 1/2 Dirac spinor for the baryon, is the Rarita-Schwinger spinor 
for the spin 3/2 baryon, v = p/M for the heavy baryon of mass M, k is the 
4-momentum of the gluon and is the corresponding propagator in axial 
gauge, r is the production vertex for the heavy quark-antiquark pair. Each 
amplitude should be 
multiplied by the color 
factor 4/3v^. 

Considerable sim- 
plification of these am- 
plitudes follows. There 
are 3 cases to consider — 
3 final state baryons. 
For the two states re- 
sulting from the vec- 
tor diquark, the elec- 
tric and magnetic am- 
plitudes must be added 
together. Then for each 
baryon, the amplitude 
is squared and a trace is 
taken to sum over spins 




Figure 2:fi{z,Q^) for 1/2' (upper pair), 3/2 (middle) 
1/2 baryons, each at Q = fi and 5.5 GeV. 
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(including spin projec- 
tion operators for the spin dependent cases) . The analog of Eq. |^ is obtained 
for each baryon. By carefully organizing the terms in the integrand, the width 
for the inclusive production of the virtual heavy quark can be divided out to 
yield the analog of Eq. ^ for each baryon. Finally the integration over s ^ 
can be performed numerically — the form factors make it difficult to write an 
analytic expression for each case. The resulting z dependent fragmentation 
functions are the "boundary" functions, obtained at a scale //^ at the thresh- 
old 4m|). To consider higher momentum scales, the Altarelli-Parisi evolution 
equations are used. 

We have taken some particular cases to illustrate the results. For the 
c(su) or c(sd) baryons, the 2c states, the diquark is given a mass of 0.95 
GeV/c^ and the ratio of diquark to hadron mass is r = 0.33. The resulting 
function, Q^) is shown in Fig. 2 for the boundary value at /x = 2md and 
for Q = 5.5 GeV. The three s-states lead to different behavior and overall 
probability. It is particularly noteworthy that the 1/2 ground state is produced 
far less frequently than the 3/2 state or the 1/2' state from the vector diquark. 
As we will see, the observed 1/2 ground states are primarily from the decays 
of these vector diquark states. 

In Fig. 3 the cor- 
responding longitudinal °" 
fragmentation function 
gi{z,Q ) is plotted for 
the two spin 1/2 states. om 
Note that it is very 
similar in shape to 
the spin averaged case. 0.02 
Lastly, the behavior of 
the transversity func- 
tion hi{z,Q'^) remains „„, 
to be determined. 

The spin depen- 
dent fragmentation func- 
tions for the spin 3/2 ' 
baryons are even richer Figure 8:51(2, Q^) for 1/2' (upper pair) and 1/2 baryons, 
in complexity. There ^^^^ at Q = fi, 5.5 GeV. 
are seven such func- 
tions at leading twist, many of which will be accessible from the decay dis- 
tributions of these states into the 1/2 state plus a pion. 
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3 Comparison with some data and summary 

The production of charmed baryons ia becoming sizeable at CESR and data 
now exist from the CLEO coUaboratioin for some of the Sc states. In particular, 
spin independent fragmentation functions have been determined for the lowest 
mass spin 1/2 and 3/2 states. The lowest 1/2+ states are the c + [u,s] and 
c+ [d, s] states, S+ and Sj^, involving the antisymmetric, spin diquarks. 
The spin 3/2+ states are c+ {u,s} and c + {d,s} baryons, S*+ and S*'', 
involving the symmetric, spin 1 diquarks. The spin 1/2+ partners, S^, of 3/2+ 
states have not been seen yet. They are presumed to have a mass below the 
Sc + TT threshold, so must be seen in radiative decay channels. Note that these 
latter 1/2+ states have the same isospin as the lower lying ground states 
1/2+ Sc and could mix with them, in principle. In any case, the measured 
fragmentation functions provide a crude test of the model. The data are fit by 
the experimenters with a common parameterization of the Peterson functioicl. 
It is easy to see in Fig. 2 that the data fall nicely on /i of our model(with 
arbitrary normalization), evolved to Q = 5.5 GeV, but are not sufficiently 
accurate to be a crucial test of the model. Note that the experimental variable 
XpU does not correspond exactly to our z, the light cone variable. 

The ratio of the |3/2 to 1/2 production can be extracted from the data 
with some uncertainty. The percentage of all 2+ states that arose from decays 
S*" — * S+ + TT^ is given as (27 ± 8)% and the percentage of all S° states that 
arose from decays S*+ ^ 2° + 7r+ is given as (17 ± 6)%. (Note that we have 
combined the statisitical and systematic errors here.) 

The experimenters do not see the 7r° channels, 2*+ 2+ + tt" and 2*° — > 
2" + tt". From isospin conservation these channels account for 1/3 of the 
decays into 2c + tt, while the reported charged tt channels constitute 2/3. 
Suppose iV2* states of both charges are produced. Then 2/3 N will be seen 
in the charged tt decay mode. The total number of 2+^°'s seen will be iV+.o — 
|Af/(0.27, 0.17) (supressing errors until the end). The number of 2+'°'s not 
coming from the decays of the 3/2 states will be N^^q — N. Assume that n+^o of 
the 2+''''s come from other fragmented states' decays. Then N+fi — N — n+fi is 
the number of direct fragmentation products of the charmed quark. The ratio 
i?(+orO) of directly fragmented 2+'° to 2*+'° is given thereby as R{+) = 
1.5 ± 0.7 - n+/N : 1 and i?(0) = 2.9 ± 1.4 - no/N : 1. 

The numbers n+.o will come from the radiative decays of the heavier 1/2 
states, as well as higher 2c states (radial and orbital excitations of the c+ {su) 
and c+ (sd) systems). We have calculated the fragmentation functions for the 
spin 1/2 quark-vector-diquark states and hence the number of 2c spin 1/2' 
states vs. 2* spin 3/2 states. That is 1.7:1 for the parameterization used in 
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Fig. 2. Assuming n+_o is due entirely to these 1/2' states decaying 100% into 
the ground state S+'°, we have for the different charge states R{+) = — 0.2±0.7 
and R{0) = 1.2 ± 1.4, both of which are consistent with the smaU ratio of 1/9 
predicted by the same model calculation. 

Hence an interesting feature of our model is that the directly fragmented 
gound state baryon is very rare compared to the vector diquark sta±^. In the 
calculations for heavy- heavy baryons by Marteynenko and SaleevEj, a ratio 
closer to unity was obtained. It will be interesting to see whether the small 
ratio, as suggested by our model, persists as more data are obtained. 
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